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Given recent indications of additional neutrino species and cosmologically significant neutrino 
masses, we analyze their signatures in the weak lensing shear power spectrum. We find that a 
shear deficit in the 20-40% range or excess in the 20-60% range cannot be explained by variations 
in parameters of the flat ACDM model that are allowed by current observations of the expansion 
history from Type la supernovae, baryon acoustic oscillations, and local measures of the Hubble 
constant //q , coupled with observations of the cosmic microwave background from WMAP9 and the 
SPT 2500 square degree survey. Hence such a shear deficit or excess would indicate large masses or 
extra species, respectively, and we find this to be independent of the flatness assumption. We also 
discuss the robustness of these predictions to cosmic acceleration physics and the means by which 
shear degeneracies in joint variation of mass and species can be broken. 



I. INTRODUCTION 

As our cosmological observations become ever more 
precise, our ability to probe smaller scales continues to 
advance, thereby allowing us to study the physics of 
structure formation beyond the standard cold dark mat- 
ter paradigm. In particular, we are now able to use 
cosmology to learn about neutrino properties, including 
the sum of their masses M^, and the effective number of 
species iVcff, both of which imprint their signatures on 
the small-scale matter power spectrum. Massive neu- 
trinos act as hot or warm dark matter, thereby sup- 
pressing structure formation below their thermal free- 
streaming scale, while adding or subtracting relativis- 
tic species changes the ratio of the acoustic and damp- 
ing angular scales of the cosmic microwave background 
(CMB) jl] [2]. The former is currently constrained to 
be > 0.05 eV by solar, atmospheric, and laboratory 
experiments [SHS], and (roughly) < 0.6 eV from cos- 
mology IBHS]. 

As for iVeff, recent oscillation and reactor experiments 
[TUl ITT] lend support to the sterile neutrino interpreta- 
tion of the LSND electron antineutrino appearance re- 
sult [12 and My > 0.4 eV while other neutrino results 
inhibit a simple global explanation (see e.g. Ref. [T3| for 
a recent review). Meanwhile, recent observations of the 
CMB damping tail (TTHH], Sunyaev Zel'dovich-selected 
clusters [IS', and baryon acoustic oscillations (BAO) [IH] 
provide further hints of extra neutrino species. 

In this paper we explore how weak gravitational lens- 
ing fits into this picture, in light of forthcoming lensing- 
optimized large-area surveys such as with the ground- 
based Dark Energy Survey (DES) [5^ and Large Synop- 
tic Survey Telescope (LSST) [2T], from the balloon-borne 
High Altitude Lensing Observatory [22], or from space 
with Euclid [53] and the Wide-Field Infrared Survey Tele- 
scope (WFIRST) [21]. Weak lensing, whereby the images 
of distant galaxies are distorted by the gravitational field 
of matter in the foreground, can be a powerful cosmolog- 
ical probe provided that we have sufficient systematics 



control. By extracting the weak lensing shear and its 
evolution with redshift we are able to robustly map out 
the gravitational potential of the Universe and how it 
changes with time. The power spectrum of this "cosmic 
shear" is directly related to the underlying matter power 
spectrum. Despite promising results, e.g. [55H25] (see 
Refs. [ 30H32j for reviews), the current constraining power 
of cosmic shear is very limited. Since we do not know the 
intrinsic shapes of individual galaxies, we must average 
them over finite patches of sky, making weak lensing a 
necessarily statistical measure whose constraining power 
is directly related to sky coverage [33l [S^ . Future data 
sets will be optimized in this respect, but in the meantime 
it is particularly timely to determine our expectations. 

We can robustly test any given cosmological model 
class by exploiting consistency relations between observ- 
ables pertaining to the expansion history and those per- 
taining to structure growth [35H38] . Given one, coupled 
with a class of cosmological models with tunable param- 
eters, we can predict the other and then compare our 
predictions to data. If the data points lie significantly 
outside of the prediction contours, then the model class 
in question is falsified. For instance, just one cluster 
that is massive and at high-enough redshift could fal- 
sify all ACDM and quintessence models if its mass and 
redshift fall significantly outside of what we predict based 
on Type la supernovae (SNe), BAO, local measurement 
of the Hubble constant (Hq), and the CMB [37] ■ In this 
way we can take advantage of the wealth of data already 
in hand from the CMB and distance measures to predict 
what we expect for these anticipated future weak lens- 
ing observations. Our analysis here builds upon [3S; to 
explore the effects of neutrinos. 

In what follows we will add neutrinos to this predic- 
tion framework, to explore how their masses and num- 
ber of species change weak lensing observables. We find 
that, for a fixed CMB, these two properties shift the cos- 
mic shear power spectrum in opposite directions; adding 
relativistic species amplifies the shear power, whereas en- 
dowing the neutrinos with nonzero masses reduces power. 
Given that quintessence can only decrease the amount of 
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power f38| , the former provides qualitatively distinct pre- 
dictions. 

This paper is organized as follows. We review our 
methodology in j jlH including the data sets we use, our 
Markov-Chain Monte Carlo (MCMC) analysis, and the 
calculation of posterior probability distributions for cos- 
mic shear observables. We then discuss our results in 
{ III including predictions for four cases - three massless 



neutrinos, three massive neutrinos, a variable number of 
massless neutrinos, and both variable mass and number. 
We conclude in fJIVl 



II. METHODOLOGY 

We describe here the data sets we use and our proce- 
dure for predicting the cosmic shear power spectrum un- 
der various assumptions for the total neutrino mass M^, 
and effective number of species A^cfF- Our methodology 
is similar to that of Refs. [35rt38j . 



A. Data sets 

We use the redshifts, luminosity distances, and sys- 
tematic uncertainty estimates of the Union2 Type la SN 
sample [39]. This sample includes 557 SNe out to a red- 
shift z = 1.12, where all light curves have been uniformly 
reanalyzed using the S ALT2 fitter [IHI ■ 

We use measurements of the BAO feature from 
Ref. [41 (which includes data from SDSS and the 2- 
degree Field Galaxy Redshift Survey), the WiggleZ Dark 
Energy Survey [42], and the SDSS-III Baryon Oscillation 
Spectroscopic Survey (BOSS) [l9l|43]. These measure- 
ments extend out to z = 0.73 and are reported as dis- 
tances relative to the sound horizon, Dv{z) /vs, where 
Dv{z) = [(1 + zfD\{z)cz/H{z)Y/^, Da is the angular 
diameter distance, H{z) is the Hubble expansion rate, 
and Ts is the sound horizon at last scattering. Even 
though these data sets have some overlap in both area 
and redshift, we treat all three as independent due to the 
different bias and type of galaxies that are targeted in 
each sample. 

Unlike Refs. [35H38] . we use CMB results from the 
most recent, 9-year release from the WMAP satellite 
(WMAP9) [H], computing the CMB angular power spec- 
tra using the code CAMB [13111]. We now further add 
the publicly available 2500 square degree release of the 
SPT measurement of the CMB damping tail [H] over the 
multipole range 650 < / < 3000, as these smaller-scale 
peaks are sensitive to neutrino physics. Per Ref. |16j . 
we treat the SZ and point-source contributions as (three) 
additional nuisance parameters, choosing the same Gaus- 
sian priors for each. 

Finally, we use the combined i?o estimate from 
Ref. E?T, of Hq = 73.8 ± 2.4 km/s/Mpc. This measure- 
ment strengthens our constraint on iVoff. 



B. Parameter sets 

We use this data set to probe flat ACDM models with 
varying fi, and iVcff. The former is the fraction of the 
dark matter density in the form of massive neutrinos 



where 



DM 



93.14/i2 eV 



(1) 



(2) 



and Ml, is the sum of the neutrino masses. The pa- 
rameter A'off is the so-called effective number of neutrino 
species gS]: 
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where pR is the energy density in relativistic species and 
pj is the energy density of photons. We increase (de- 
crease) A^eff by adding (subtracting) massless species. 
The default is three species with degenerate masses. Un- 
fortunately this means that A^eff < 3.046 leads to a nega- 
tive number of massless species, an unphysical situation 
which is treated by CAMB as a negative energy den- 
sity; we follow the SPT analyses [T5HI8] and ignore this 
since A^^ff > 3.046 is highly favored and indeed these un- 
physical cases make up no more than a few percent of 
the samples in our MCMC chains. We then set the pri- 
mordial helium abundance Yp from the physical baryon 
density fl\^h^ by requiring "big bang nucleosynthesis con- 
sistency" [in], such that 

Yp = 0.2485 -I- 0.0016 [273.9f7b/i^ - 6 + 100 (5 - 1)] , 



where S depends on the number of neutrino species 

5^ = 1 + (A^cff - 3.046) . 
The total parameter set we use is 



(4) 



(5) 



6> = {rJb/i^ nuuh^, T, 0A, n,, In^, TVeff} , (6) 

where flumh'^ is the present physical dark matter den- 
sity relative to the critical density, r is the reioniza- 
tion optical depth, 9a is the angular size of the acous- 
tic scale at last scattering, is the spectral index of 
the power spectrum of initial fluctuations, and As is 
the amplitude of the initial curvature power spectrum 
at fcp = 0.05 Mpc~^. All other parameters, including 
the Hubble constant — 100 ft. km/s/Mpc, the present 
total matter density flm, the dark energy density fluE, 
and the amplitude of the matter power spectrum today 
as, can be derived from this set. We will study four dif- 
ferent cases: (1) fu = and A^off — 3.046, i.e. standard 
flat ACDM, (2) is allowed to vary A^^ff = 3.046, (3) 
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Ncs is allowed to vary, fi, = 0, and (4) both and A'cff 
allowed to vary. 

For a given set of parameters 9 that defines the cos- 
mological model class in question, we use the CosmoMC 
code [SOI [ST] to sample from the joint posterior distribu- 
tion, 



P(0|x) = 



J d9C{x\9)V{e) 



(7) 



where £(x|0) is the likelihood of the dataset x given the 
model parameters 9 and V{9) is the prior probability 
density. For the standard ACDM parameters we use the 
same priors as in |38| (flat priors that are wide enough 
to not limit our constraints), and we similarly choose 
wide priors for f^, and A^eff that are informed by the cur- 
rent limits from data as summarized in Refs. [HI l52] . 
In particular we choose the extremely conservative prior 
1.047 < iVeff < 10.0, as in the SPT analyses. 



C. Weak lensing observables 

We can compute the posterior probability distribution 
for any derived statistic from the joint posterior distri- 
bution of the cosmological parameters. In particular, in 
order to compute the cosmic shear power spectrum, we 
must first compute the comoving angular diameter dis- 
tance D and the nonlinear matter power spectrum A^^. 
In a flat universe (curvature f^K = 0) , the former is equal 
to the comoving radial coordinate and is related to the 
cosmological parameters through 



D(z) 



dz' 



Here the Hubble expansion rate is 

H{z) = Ho [a„(i + + (1 - a^)] 

where the total matter density is 



DM 



(8) 



(9) 



(10) 



and the contribution from radiation is assumed to be 
negligible. 

We compute the z = linear matter power spectrum 
AL(fc;0) using CAMB. The linear matter power spec- 
trum at earlier redshifts then depends on the growth 
function of linear density perturbations. Massive neu- 
trinos suppress growth in a scale-dependent manner, and 
we model this using the Eisenstein and Hu [S3] fitting 
function. 



Alik;z)^Al{k;0) 



T\k,z) Dljz) 



(11) 



where T{k,z) and Di{z) are given by their Eqs. (7) and 
(8), respectively, which we have modified according to 



Ref. [51] to improve accuracy in the case of three mas- 
sive neutrinos. Note that Di{z) corresponds to the stan- 
dard scale- independent growth function in the absence of 
neutrinos, and the scale-dependent effects of their free- 
streaming are encoded in T{k,z). Comparing to results 
from CAMB for nonzero redshifts, this fitting formula 
typically reproduces the growth to better than 1% for all 
k and z we use here to compute the shear. 

We compute the full nonlinear matter power spectrum 
at a given redshift using the Halofit fitting function [SS] 
(see Ref. [35] for a summary), modified for the effects of 
massive neutrinos [52) . The original Halofit fitting func- 
tions have been found to only be accurate (even for the 
flat ACDM model) at up to the 5-10% level compared 
with A^-body results, for instance with the Coyote Uni- 
verse project ^56-r.-58j . We flnd that whether or not we use 
the massive neutrino modification [52] leads to errors of 
order a few percent for the small neutrino masses consid- 
ered here. These systematic errors are smaller than the 
statistical errors arising from our data sets in the same 
regime [38], and so we expect our results to be fairly 
robust to them. Likewise, for a wide range of baryonic 
effects, systematic shifts are at most comparable to cur- 
rent statistical errors [38^. 

The shear (or equivalently the convergence) power 
spectrum is then equal to 



2^ -W^'^^ 



dz^mA^ 



H 



NL 



^;-) > (12) 



where k « l/D in units of Mpc~^ in the Limber ap- 
proximation and we have defined the geometric lensing 
efficiency factor 



dz'n{z') 



D' - D 
D' 



(13) 



The efficiency factor weights according to the source dis- 
tribution in a given survey, n{z), normalized such that 
n{z)dz = 1. Here we use the model 



n(z) (X — 
\2o 



exp 



(14) 



with parameters (zo,a,/3) — (0.555,1.197,1.193) for a 
simplified model ground-based survey, such as CFHTLS 
or DES, with an approximate median redshift of 0.8. 



III. RESULTS 



The histograms in Fig. [T] illustrate what happens 
to our predictions for our six flat ACDM parameters 
{nbh'^,nuMh'^,T,0A,ns,lnAs,f^,Nes} when and/or 
A'off are allowed to vary. All of our results from here on 
are presented with the following color-coding: 



(1) Blue: 
ACDM 



and A^, 



cff 



3.046, i.e. standard flat 
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FIG. 1. One-dimensional constraints on the six cosmological parameters of our baseline fiat ACDM model for our four cases: (1) 
three massless neutrinos (blue, thick solid); (2) f^, allowed to vary (magenta, dashed); (3) A^cft allowed to vary (red, dot-dashed); 
and (4) both and A'cff allowed to vary (green, thick dotted). 



(2) Magenta: is allowed to vary and N^ff = 3.046 

(3) Red: A'cff is allowed to vary and fi, = 

(4) Green: both and A'cff are allowed to vary 

We further summarize our constraints on and A'cff 
in Table |T] for each of these four cases. In the most gen- 
eral case (4), we constrain M^, < 0.67 eV at the 95% 
confidence level and we find A'cff = 3.71 ± 0.35. These 
constraints are in agreement with the current state-of- 
the-art as seen in the literature for cosmological probes, 
e.g. Ref. [59j. On the other hand, since the M^, con- 
straints are based on the CMB and expansion history 
measurements rather than growth measurements (see e.g. 
[60]), they are less robust to generalizations of the flat 
ACDM model, e.g. the addition of spatial curvature (see 
below). 

TABLE I. Constraints on the sum of the neutrino masses 
Mj^ and effective number of species A'cff for our four cases: 
(1) three massless neutrinos; (2) allowed to vary; (3) A'efr 
allowed to vary; and (4) both and A'cff allowed to vary. We 
report the 95% upper limit on M^, and the mean and 68% 
confidence interval about the mean for A'cff. 





1 


2 


3 


4 


(eV) 




< 0.45 




< 0.67 


A'cff 






3.56 ±0.32 


3.71 ±0.35 



We find that allowing these small neutrino masses does 
not significantly change any of our parameter constraints, 
whereas allowing additional sterile neutrino species does. 
In particular, we see the physical dark matter density 
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FIG. 2. Constraints in the SIdm/i^ — A^eff (left panel) and 
ns — TVcff (right panel) planes without (yellow) and with (grey) 
SPT CMB data [TS], for the case where both A'cff and fi, are 
varied. 

^^DM^^ is increased while the power spectrum of primor- 
dial fiuctuations gains an enhancement in both the tilt 
Us and amplitude Ag [61] . 

In Fig. [2] we show how the degeneracies between A'cff 
and r^DM" or tighten with the addition of the 
SPT CMB data, which provides several more peaks in 
the small-scale regime. As has been noted elsewhere, 
e.g. |59| . we find that the inclusion of an Hq prior 
strengthens our constraints on A'cff. However, we still get 
meaningful results without it due to our inclusion of SPT 
data. The damping scale test provides constraints which 
are independent of low-redshift dynamics, and therefore 
the specifics of the dark energy model. 

In Fig. [3] we show how the z = linear matter power 
spectrum prediction contours shift with the addition of 
massive neutrinos, additional neutrino species, or both. 
The top-left panel shows the baseline (i.e. three mass- 
less neutrinos) prediction, and the top-right and bottom 
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FIG. 3. Flat ACDM predictions for the z — Q linear matter power spectrum, with k in units of h Mpc~^ and P{k) in units 
of Mpc^, and the color-coding as before for the four cases: (1) three massless neutrinos (blue); (2) allowed to vary 
(magenta); (3) A^'cff allowed to vary (red); and (4) both and A'cff allowed to vary (green). The top-right and bottom panels 
are all plotted with respect to the maximum likelihood "blue" model prediction, showing the 68% and 95% confidence level 
regions, and with the same axis scales for comparison. 



panels are plotted with respect to the maximum like- 
lihood baseline model prediction, with color-coding as 
before. For plotting purposes we follow the usual con- 
vention of taking P{k) — {2tt^ / k^) (k) , with k in units 
of h Mpc~^. We see that endowing neutrinos with mass 
serves to suppress structure growth, in accordance with 
conventional wisdom, despite the similar parameter pre- 
dictions as seen in Fig. [T] We further see that allowing 
for additional sterile neutrino species serves to enhance 
structure on small scales. This is because of an N^g — Ug 
degeneracy. A larger A^off suppresses power in the high £ 
CMB spectrum due to damping which then allows a com- 
pensating increase in rig or the high-fc primordial power 
spectrum. Given the preference for additional species 
seen in Table |lj we find that this significantly shifts our 
Pl contours up for large k. 

We show the resulting 2D cosmic shear power spectra 
in Fig. |4| for our model ground-based weak lensing sur- 
vey, again with the top-left panel showing the baseline 
prediction, and the top-right and bottom panels plotted 
with respect to the maximum likelihood baseline model 
prediction. 



In the context of the flat ACDM model, an observed 
deficit of small-scale cosmic shear of between 20 — 40% 
would indicate finite neutrino mass and could not be ex- 
plained by other currently allowed cosmological param- 
eter variations. An observed excess of 20 — 60% would 
indicate extra neutrino species and comes from the free- 
dom to raise the tilt due to the A'cff — rig degeneracy in 
the CMB. Cosmic shear measurements provide a means 
of breaking this degeneracy. We illustrate the issue in 
Fig. |5j where we plot our constraints on the cosmic shear 
power at I — 1000 vs Ug for the case where A'cff is var- 
ied but fi, = 0. The case where both f^, and A'cff are 
allowed to vary is harder to distinguish in that the two 
effects can partially compensate for each other. On the 
other hand, a further breaking of the A^efi ~ ng degen- 
eracy is expected from the Planck survey |62j . thereby 
allowing these mixed cases to be better separated with 
cosmic shear. 

It is interesting to note that while the addition of SPT 
data tightens the predictions for the shear power spec- 
trum it actually weakens and shifts the predictions on 
the growth function, as we show in Fig. [6] for 2 = 0. 
This reflects a mild tension between this data set and the 



6 



Vary only /v 




10 100 1000 10 100 1000 

/ / 



FIG. 4. Flat ACDM predictions for the ground-based cosmic shear power spectrum, with the color-coding as before for the 
four cases: (1) three massless neutrinos (blue); (2) f^, allowed to vary (magenta); (3) A'cff allowed to vary (red); and (4) both 
and A'eff allowed to vary (green). The top right and bottom panels are all plotted with respect to the maximum likelihood 
"blue" model prediction, showing the 68% and 95% confidence level regions, and with the same axis scales for comparison. 
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0.94 0.96 0.98 1 

FIG. 5. Constraints on the ground-based cosmic shear power 
spectrum amplitude at / — 1000 (multiplied by 10^") vs Us, 
for the case where A^cff is varied but — 0, showing the 68% 
and 95% contours. 



BAO measurements caused by its improved measurement 
of f^Divi/i^- Unlike the similar tension between BAO and 



Hq for the flat ACDM model, this tension is not allevi- 
ated by allowing A'eff (or /^) to vary. On the other hand, 
constraints on the growth function are not the dominant 
source of error for shear predictions and so this tension 
is not relevant for our purposes. 

Beyond the flat ACDM model, there are other possi- 
bilities that can explain a deficit or excess of small-scale 
shear. When the dark energy equation of state is general- 
ized to allow quintessence, only a deficit can arise due to 
the restriction that w > —1 iMj- Hence these cases can 
masquerade as massive neutrino models unless further 
information on the shape and redshift dependence of the 
power spectrum is obtained. Without neutrino number 
changes, an excess cannot be explained by quintessence 
and hence would indicate more exotic cosmic acceleration 
physics with enhanced forces in the dark sector. Again, 
the expected improvements from the Planck survey will 
help distinguish between these possibilities. 

Likewise we have also tested the robustness of these re- 
sults to dropping the flatness assumption. Allowing for 
curvature significantly degrades our constraints on M^, 
where we find the 95% limits expand to M^, < 1.68 eV 
and < 1.88 eV for the cases when f^, only is varied 
and for when both and A^off are varied, respectively. 
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^1^=0) 

FIG. 6. One-dimensional constraints on the scale-independent 
growth function D\ aX z = 0, for the case where both A^'cff 
and fv are allowed to vary, without (black solid) and with 
(green dotted) the addition of SPT CMB data. Note that 
the additional data shifts and weakens the joint constraint, 
thereby indicating tension in the data sets. 

On the other hand, our constraints on A'cfr are not signif- 
icantly different then those of the flat case. We further 
find that our high-redshift growth function constraints 
are weakened, but since the confidence contours are still 
well within the 1% range we find that there is not a sig- 
nificant effect on the shear predictions. Indeed our shear 
predictions are qualitatively the same. 

IV. DISCUSSION 

Inspired by recent evidence for massive neutrinos and 
the possibility of additional species, we have provided an 
analysis of the signatures of such "nonstandard" neutrino 
physics on the weak lensing shear power spectrum. By 
using observations of the expansion history from Type la 
SNe, BAO, and local measures of i7o, coupled with ob- 
servations of the CMB, we can predict future structure- 
growth observables such as those from weak lensing. 
From doing so for our four different scenarios - the stan- 
dard case with three massless neutrinos, three massive 
neutrinos, any number of massless neutrinos, and three 
massive neutrinos with any number of massless neutrinos 
- we can look for signatures that cannot be mimicked by 
any currently allowed variation in the other parameters 
of the flat ACDM model. We present only results for a 



representative model ground-based survey here, but the 
results for any survey configuration will be qualitatively 
the same. 

For our most general case, where we vary both neu- 
trino mass and number of species, we find parameter 
constraints that are consistent with the current litera- 
ture. Using only distance measures and the CMB, we 
are able to constrain M^, < 0.67 eV at the 95% confi- 
dence level and N^ff — 3.71 ± 0.35. 

Such variations from the standard neutrino parameters 
allow changes in the predicted shear power spectrum that 
cannot be mimicked by other flat ACDM parameters. For 
example shear deficit in the 20-40% range would indicate 
neutrino masses near saturation of current bounds with 
A'cfi ~ 3 whereas shear excess in the 20-60% range would 
indicate extra neutrino species and lower masses. The 
latter is because of a partial degeneracy between raising 
iVoff and spectral tilt of the primordial power spectrum 
Us (see Fig. [5]). 

It has been noted previously that generalizing from the 
flat ACDM model to quintessence can only serve to re- 
duce the matter power spectrum 37, 38^ , within the con- 
text of the standard (three massless) neutrino scenario. 
Quintessence effects can therefore mimic the deficit pre- 
dicted by massive neutrinos without further information 
from the shape and redshift dependence of the power 
spectrum. 

A future measurement of excess cosmic shear could 
provide supporting evidence for extra neutrino species. 
An excess could also be explained by more exotic acceler- 
ation physics that enhances structure growth. Since the 
neutrino effect comes from the primordial power spec- 
trum, improved measurements from the Planck satellite 
would help distinguish these options should an excess be 
found. 
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